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ABSTRACT 

We study the observed evolution of galaxy clustering as a function of redshift. We 
find that the clustering of galaxies, parameterized by the amplitude of fluctuations 
in the distribution of galaxies at a co- moving scale of Sh"^ Mpc, decreases as we go 
from observations of the local Universe to z ~ 2. On the other hand, clustering of the 
Lyman break galaxies at 2; ~ 3 is very strong, comparable to the clustering of present 
day galaxies. 

However there are three major factors to take into account while comparing clus- 
tering measurements coming from various surveys: the so-called "scale-dependence" 
effect, due to measurements being made at different scales; the "type-selection" effect 
introduced by the fact that different galaxy surveys select different populations which 
do not have the same clustering amplitudes; and the Malmquist bias which means that 
within a given survey the more distant galaxies tend to have brighter absolute magni- 
tudes, and so do not have the same clustering amplitude. We correct for the first two 
effects and discuss the implications of Malmquist bias on the interpretation of the data 
at different z. Then we compare the observed galaxy clustering with models for the 
evolution of clustering in some fixed cosmologies. Correcting for the scale-dependence 
effect significantly reduces the discrepancies amongst different measurements. 

We interpret the observed clustering signal at high redshift as coming from objects 
which are highly biased with respect to the underlying distribution of mass; this is not 
the case for z < 2 where measurements are compatible with the assumption of a much 
lower biasing level which only shows a weak dependence on z. Present observations 
still do not provide a strong constraint because of the large uncertainties but clear 
distinctions will be possible when larger datasets from surveys in progress become 
available. Finally we propose a model-independent test that can be used to place a 
lower limit on the density parameter fio- 
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1 INTRODUCTION 



hav 1998; Tegmark and Peebles 1998; [Narayanan, Berlind 



and clusLeis ul 



It IS believeu Lliat stiacLLiies liJ^e galaxies and clusteis 
galaxies formed by accretion of matter onto small inhomo- 
geneities present in the early Universe. The simplest models 
assume that the distribution of galaxies is directly related to 
the underlying density distribution and the two distributions 
evolve in a similar manner. This has provided a key motiva- 
tion for redshift surveys of galaxies. However, many studies 



and Weinberg 1998| ; [Blanton et al. 19981 Folin et al. 199^ 



have shpwii UiaL Lliti relation bfelweeii galaxy cluHlei'liig and 



that of midetlyiiig maiiei is not siiiiple and this lelatioii is 



in genera l, a function of time (Brainerd and Villuinscn 1994t 
Fry 1996|; |Mo and White 1996t [Bagia 1998|; [Dekel and La 



Baugh et al. 1999; Kauffmann et al. 199!:). These studies 



deal with the evolution of the clustering of dark matter ha- 
los and many factors, including observational selection func- 
tions and evolution of stellar populations in galaxies have to 
be taken into account before these results can be applied to 
real observations. Some comparisons of models and observa- 



tions have been carried out (Matarrese et al. 1997; Moscar- _ 
dini et al. 1998|) where a large class of models were compared 



with the available observations. 



Models of halo clustering (Mo and White 1996; Bagla 
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1998) suggest that more massive halos cluster more strongly 
than halos o f lower mass . This follows from the simple model 
for biasing (Kaiser 1984) in which rare objects cluster more 
strongly than the more typical objects. If, as observations 
suggest, the luminosity of a galaxy increases monotonically 
with the mass of the halo in which it resides, then we expect 
brighter galaxies to cluster more strongly than fainter galax- 
ies. The same argument also suggests that amongst halos of 
a given mass, older halos cluster more strongly. If the grav- 
itational clock is synchronised with the stellar clock then 
we expect early type galaxies to cluster more strongly than 



late typ e galaxies. This is seen in simulations ( Blanton et al 



199S) t lat include simple recipes for star formation. 



A common conclusion of all theoretical and numeri- 
cal studies of halo clustering is that the rate of evolution 
of halo clustering, Dhaioii), is always slower than the rate 
of evolution of clustering in dark matter Dm{i), so that 
Dhaio{t) < Dm{t), where the dot represents differentiation 
with respect to time. These rates are equal only in the limit 
when all the matter has collapsed into halos. 

Turning now to a brief discussion of effects that influ- 
ence the observed evolution of galaxy clustering, consider 
a Universe in which the galaxies do not evolve: neither in 
their stellar content, nor in their distribution in space. So 
the clustering is fixed in comoving space and galaxies at all 
redshifts are similar to the ones we see in the local Universe. 
What will be the observed amplitude of clustering at dif- 
ferent redshifts, if we conduct a magnitude limited redshift 
survey in such a Universe? As nothing is changing as far as 
galaxies are concerned, the only differences are given by ob- 
servational selection effects. There are at least two of these: 
Malmquist bias and K-correction. Let us examine the effect 
of these two factors separately. 

In an apparent magnitude limited survey, we will only 
see brighter galaxies at high redshifts whereas at lower red- 
shifts we will also see fainter galaxies. Since brighter j 



tend to cluster more strongl y than fainter galaxies (Park et 
al. 1994; Loveday et al. 1995| ) , the effect of Malmquist bias in 
our imaginary survey will lead to an apparent increase in the 
amplitude of clustering with redshift. Also in an apparent 
magnitude limited sample, the observed clustering ampli- 
tude will always exceed the true clustering amplitude of all 
galaxies at that redshift. This variation with redshift can 
be avoided by using a cut in absolute magnitude instead of 
apparent magnitude, but the inferred absolute magnitudes 
depend explicitly on the assumed values of cosmological pa- 
rameters. However, without going into any details of depen- 
dence on cosmological parameters and models of galaxy evo- 
lution, we can conclude that the observed rate of evolution 
of clustering amplitude is always smaller than the true rate 
in a magnitude limited survey, the only assumption here be- 
ing that brighter galaxies cluster more strongly than fainter 
galaxies: boba{t) < Dtrue{t). 

The effect of K-correction depends strongly on the 
wave-band used to define the sample. For example, the dif- 
ference between the rest-frame B band luminosity for ellipti- 
cals and irregulars observed in the B band is ~ 4 magnitudes 
at z = 1, but is only ~ 2 magnitudes if they are observed 
in the I band (Tresse 1999). This effect means that the rela- 
tive populations of early- and late-type galaxies in a B band 
selected survey will vary strongly as a function of redshift, 
preferentially selecting late-types at higher redshift. Since 



later-type galaxies cluster less strongly than early types (e.g. 
Loveday et al. 1995; Hermit et al. 1996; Guzzo et al. 1997; 
Loveday Tresse & Maddox 1999) a sample defined in the 
optical will generally tend to underestimate the correlation 
amplitude at higher redshifts. In the near infra-red the spec- 
tral energy distributions of different galaxy types are very 
similar to each other, and so the K-corrections are very simi- 
lar. Hence the problem is less significant if we use near infra- 
red wavelengths to define the sample. Quantifying this ef- 
fect requires knowledge of the relative fractions of different 
galaxy populations as a function of redshift, and of how the 
clustering varies for each population. We plan to tackle this 
problem in a future paper, but in our present analysis we 
have simply grouped the various measurements according to 
the survey selection criteria so that the differences between 
types are minimized. 

This problem could be by-passed if a fixed rest-frame 
bandpass were used to define the sample, but even this would 
not take into account the spectral evolution of galaxies. In 
the real Universe, all galaxies tend to get bluer at higher red- 
shifts, because of stellar evolution, and so the differences be- 
tween different types of galaxies become smaller. This means 
that in the real Universe, uncertainties due to different K- 
corrections will be smaller than in our imaginary Universe. 

Another factor which affects the interpretation of the 
data from various surveys is the "scale-dependence" of the 
measurements: different surveys sample clustering on differ- 
ent physical scales. Even though the values for the correla- 
tion length ro{z) may be of the order of a few Mpc's, many 
surveys (especially those at high redshifts) do not sample 
scales above ~1 h~^ Mpc. Since the correlation function is 
not necessarily a pure power law and the measurements have 
different slopes, this effect can introduce systematic offsets 
between the various surveys. To take this into account in this 
paper, we always compare theoretical predictions to data at 
the same scale. 

This paper presents a compilation of clustering mea- 
surements and aims to introduce some simple scenarios for 
the evolution of clustering in order to understand the dif- 
ferent factors which enter into the comparison of data and 
models. More detailed quantitative discussion on the evolu- 
tion of clustering will require corrections for all of the effects 
discussed above, possibly through comparison with a library 
of models for galaxy formation (e.g. semi-analytical models) 
convolved with observational selection effects. 

The layout of the paper is as follows: section 2 intro- 
duces simple analytical models for the evolution of bias with 
epoch, while section 3 presents the observations and shows 
the trend of the "raw" data (i.e. with no corrections ap- 
plied for the effects described in this section) together with 
simple models for the evolution of clustering. A more de- 
tailed analysis of the data, together with our interpretation 
of the results is given in section 4. Section 5 introduces a 
model-independent test to place a lower-limit on the density 
parameter no, while section 6 summarises our conclusions. 



2 EVOLUTION OF BIAS 

Most theoretical models for evolution of galaxy clustering 
identify galaxies with halos. Here we will describe a few 
models for the redshift evolution of the bias of halos b{z). 
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• No Evolution (BO) : In this model bias does not 
change with redshift and remains constant at its present 
value. This assumption does not have any physical basis and 
this model serves only as a reference. The bias is defined as 



feo = 



CTS,! 



(1) 



where as is the r.m.s. density fluctuation at 8/i~^ Mpc, sub- 
scripts g and m denote galaxies and total underlying mass, 
respectively and all quantities are evaluated a.t z = Q. 
As the halo correlation function evolves at a slower rate 
than the dark matter one, we expect bias to be higher at 
high redshifts than at present. Therefore this model should 



termined by the choice of the cosmological parameters and 
power spectrum of mass fluctuations. In the following mod- 
els we will consider three combinations, namely 

• I : ^0 = 1 ^A.o =0 /lo = 1 r = 0.25 
4% = 0.65 

• II : fio = 0.4 VLk.o = fto = 0.65 F = 0.23 
0^%, = 0.64 

• III : 9,0 = 0.4 ^A.o = 0.6 ho = 0.65 F = 0.23 
4'.™ = 1-07 

The normalization of the linear r.m.s. mass density fluctua- 
tions at d,h~^ Mpc (Tg"J„ for each power spectrum i s chosen to 



underetiLimaLe iis,g{z) &i high 



match the f our year COBE DMR observations (Bunn and 



• T est Pai ' ticl i s Bias (Bl) . This model does not as- 7: -r.. r ,■ r i-rr . 1 . 1 ,-r 

...... ,... equation y3) as functions of dinerent cosmologies and du- 

mp fiTn/tnincr nnniif f np nriffm rtr VipilnQ nr nnniit f ripir ini- ^— ' 



sume anything about the origin of halos, or about their ini- 
tial distribution. Test particles are distributed through the 
Universe such that their density contrast is proportional to 
the density contrast of the total underlying mass. The model 
describes the evolution of bias for these test particles by as- 
s uming that they follow t he cosmic flow. It can be shown 



White 1997). The values of the parameters 6g and /3 in 



( Nusser and Davis 1994 Fry 1996; Tegmark and Feeble; 
199S ) t lat the bias for such test particles evolves as 



b{z) = 1 -I- 



60 



D{z 



(2) 



where D{z) is the linear growth rate for clustering. Here 60 
is the bias for the set of halos/test particles at present epoch 
(see equation hi). This bias does not depend on the mass of 
halos and the mod el works we ll in the range < 2: < 1 for 
CDM hke models ( JBagla 1998| ) if the halo distribution is bi- 
ased. Note that the predicted variation for anti-biased halos 
(&o < 1) is not seen in simulations simply because the ba- 
sic premise of inert, indestructible halos is not co rrect. This 
model is also called the galaxy conserving model (Matarrese 



et al. I!p97|). 

• Merging Model (B2) : This model (Matarrese et al 



i)/^(^f, 



(3) 



1997) f c r the evolution of galaxy bias allows halos to undergo 
dissipative collapse (i.e. merging). It is based on the Mo and 
White (1996) model for halo bias which was in turn com- 
puted, from the formalism of Press and Schechter (1974). 
The bias is computed for all halos with mass M above a 
certain threshold Mmin which is computed by normalizing 
the effective bias b"^" at z = to the observed bias. The 
generic expression for bias of halos of mass A/ > Mmin is 
given by 

b{M,z) = b-^ + {bf ~ 

where the parameters b^ and /3 depend on the choice of 
Mmin and the background cosmology. As in the model of 
Mo and White (1996), there is also some dependence on the 
formation redshift of halos, which in this case was fixed to 
a generic zp- The constant value 6_i = 0.41 is obtained in 
the limit Mmin — > 0. 

In the present study, we decided to leave Mmin as a free 
parameter (ranging 10^° M© < Mmin < 10^"^ M0) so as to 
allow the objects under consideration to be either transient 
or to h ave changed their prop erties with time (c.f. transient 
model (Matarrese et al. 1997)), and to keep the discussion 
on a more general level. 

As already stated, both the values of 6q and (3 and those 
for bo (and therefore of as,m through equation hi) are de- 
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ferent halo masses are given in Matarrese et al. (1997) and 
Moscardini et al.(1998). However, much of the discussion is 
independent of the specific normalization or the shape of the 
dark matter power spectrum. 

We note that biasing is likely to be non-linear, non-local, 
scale-dependent, type-dependent and stochastic (e.g. Dekel 
& Lahav 1998; Tegmark & Peebles 1998; Blanton et al. 1998; 
Narayanan, Berlind and Weinberg 1998), so the models dis- 
cussed here will give a highly simplified picture of galaxy 
clustering evolution. 



3 OBSERVATIONS 

In this section we put together the observations of galaxy 
clustering from different surveys. In the following we assume 
that the correlation function has a power-law form ^(r, z) = 
[r/ro[z)]~~' at the relevant scales, and hence the index 7 
and the correlation length ro are sufficient to describe it at 
a given redshift. Most observers quote these numbers and 
in cases where these were given in proper coordinates, we 
converted them to the corresponding comoving scales. To 
transform results between different cosmologies, we use the 
following expression 



'■02(2) = 



ftpi 

ho2 



Xljz) 
X2{z) 



1-7 



P{^oi,z) Fijz) 

P{no2,z) F2{z) 



1/7 



m(2).(4) 



This equation is derived by requiring that the angular cor- 
relation function of a set of galaxies between redshift z and 
z + Az is the same in different cosmologies. If Az is small 
enough for us to assume a constant redshift distrib ution of 
objec ts N{z) then the relativistic Limber equation ( Peebles 



198C ) leads us to equation H for a power law correlation func- 
tion. The relative expressions for the comoving coordinate 
X and the functions P and F for different geometries are 
given in Appendix 1 (see also Magliocchetti and Maddox, 
1999, Treyer and Lahav, 1996). 

To compare different datasets, we will use the r.m.s. 
ffuctuations in the galaxy distribution at the scale of 
8/1"^ Mpc, 



as. 



We relate as to the other two parameters as (Peebles 



198C) 



0-8(z) 



ro{z) 



-,1/2 



(5) 
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Figure 1. Measurements of eg as a function of redsliift respectively for an Einstein-deSitter Universe Co = ^, ho = 1 (left panel), OCDM 
Qo = 0.4, ho = 0.65 (middle panel) and ACDM Oq = 0.4, H^.O) ho = 0.65 (right panel). The error bars are the quoted 68% confidence 
limits or the equivalent taken from the papers cited in the description of surveys used here. The thick line shows the linear rate of growth 
for dark matter. The dashed curves show evolution as predicted by the test particle model, the two curves are anchored to the APM and 
the IRAS observations at z = 0. The dotted curves show the evolution of bias in the merging model for 10^'^, lO^'^, lO^'^ and lO^'^M 
from bottom upwards. 







72 



(3-7)(4-7)(6-7)2^' 



(6) 



The list of surveys from which the data points have been 
taken is given below and in table 1, along with detailed infor- 
mation about the relevant parameters of each survey. These 
parameters include the median redshift (z) of the survey, 
the number of galaxies Nqal in each survey, the angular 
coverage, the selection band, the range of absolute magni- 
tudes L (if available) of the objects included in the analysis 
of clustering and the relative values of 7. We considered the 
following surveys: 



• Stroml o-APM (|Loveday et al. 1995| ). 

• IRAS (ISaunders et al. 1992|). 

• Las Ca mpanas (|Huan et al. 1996 ). 

• CFRS ( [LeFevreet al. 199^ ). 

• HDFl (Connolly, Szalay & Brummer,1998). 

• HDF2 (Maglio cchetti fc Maddox, 1 999). 

• Keck K- band (jCarlberg et al. 1997| ). 

• CNOC2 dCarlberget al. 19981 ). 

• LBGi (Giavalisco et al., 1998). 

• LBG2 (Adelberger et al., 1998). 

The values for as{z) for the three cosmologies presented in 
section 2 are shown in Figure 1. Note that there are two 
different points for Lyman break galaxies (LBG) at 2; = 3. 
These correspond to two different subsamples - one with 
observed redshifts that can generally be described as the 
brighter sample, and the other with photometric redshifts. 
In the first case th e amplitude of fluctuati ons is determined 
by counts-in-cells ( Adelberger et al. 1998 ) and in the other 
case it is arrived at thr ough the angular correlation function 
(Giavalisco et al. 1998). The higher amplitude point corre- 



In order to guide the eye we also plotted some simple models 
for the evolution of clustering. Here we assume o"8,m(z) to 
vary with redshift according to linear theory, i.e. we write 



0-8,9(2;) = b{z) ag^rniz), with 



Crs,miz) = (7s^m{z = 0) 



D(z = 0)' 



(7) 



{D(z = 0) = 1 by deflnition), and b{z) given by equations 
(^j^. The thick line shows the linear rate of growth for 
dark matter. It is clear that within any given survey, the 
amplitude of fluctuations does not fall as rapidly as the lin- 
ear rate. This is encouraging because according to the ar- 
guments outlined in the introduction, the observed rate of 
evolution should always be smaller than the rate of evolution 
of mass. The dashed curves show evolution as predicted by 
the test particle model (Bl), the two curves are anchored to 
the APM and the IRAS observations at z = 0. The dotted 
curves show the evolution of bias in the merging model (B2) 
for 10", 10", 10^^ and lO^^M© from bottom upwards. 
The basic pattern followed by the amplitude of clustering, 
even though it is masked to some extent by large error bars 
and differences in different datasets, is that - independent 
of the cosmological model - at low redshifts, as,g decreases 
with increasing redshift, reaches a minima around z — 2 and 
then rises again at higher redshifts. This type of variation 
has been seen for dark matter halos in N-body simulations 
(see e.g Jenkins et al. 1998), but given the observational 
complications discussed earlier, these simulations cannot be 
directly compared to the observational data. 



sponds to the sample with redshifts. 
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Table 1 . Summary of the properties of each survey. 
Survey (z) l^GAL Angular Coverage Selection band 



Absolute Magnitude 



IRAS 


<0.1 


9080 


|b| > 10° 


Ffio < OSJy 


7 






~0.06 


1.57 ±0.03 


Stromlo-APM 


0.06 


1757 


4300 deg2 


bj < 17.5 




-22 < Lb < -15 


~0.05 


1.71 ±0.05 


Las Campanas 


0.1 


19558 


700 deg2 


R < 17.75 




-23 < Lr< -18 


0.13 


1.85 ±0.05 


CNOC2 


0.36 


~2000 


40 X 9' X 8' 


R< 24 




Lr 


> 


-20 


0.08 


1.8 ±0.1 


CNOC2 


0.36 


~2000 


40 X 9' X 8' 


_R< 24 




Lr 


> 


-20 


0.28 


1.8 ±0.1 


CNOC2 


0.36 


~2000 


40 X 9' X 8' 


ii;<24 




Lr 


> 


-20 


0.43 


1.8 ±0.2 


CFRS 


0.56 


591 


5 X 10' X 10' 


17.5 < Iab < 22.5 


-21.5 


'^Lb< -18.5 


0.35 


1.64 ±0.05 


CFRS 


0.56 


591 


5 X 10' X 10' 


17.5 < Iab < 22.5 


-22.0 


<Lb< -19.5 


0.62 


1.64 ±0.05 


CFRS 


0.56 


591 


5 X 10' X 10' 


17.5 < Iab < 22.5 


-22.5 


<Ls<-21.0 


0.87 


1.64 ±0.05 


Keck K-band 


0.7 


248 


27 arcmin-^ 


K <20 




Lk 


> 


-21.5 


0.34 


1.8 


Keck K-band 


0.7 


248 


27 arcmin^ 


K <20 




Lk 


> 


-23.5 


0.62 


1.8 


Keck K-band 


0.7 


248 


27 arcmin^ 


K <20 




Lk 


> 


-23.5 


0.97 


1.8 


Keck K-band 


0.7 


248 


27 arcmin^ 


K <20 




Lk 


> 


-23.5 


1.39 


1.8 


HDFi 


1 


926 


5 arcm 


n2 


Iab < 27 




7 






0.6 


1.8 


HDFi 


1 


926 


5 arcm 


n2 


Iab < 27 




7 






1.0 


1.8 


HDFi 


1 


926 


5 arcm 


n2 


Iab < 27 




7 






1.4 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






0.6 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






1.0 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






1.4 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






1.8 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






2.2 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






2.6 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






3.0 


1.8 


HDF2 


1.6 


946 


4 arcm 


n2 


AB(8140) 


< 28 


7 






3.6 


1.8 


HDF2 


1.6 


946 


4 arcmin^ 


AB(8140) 


< 28 


7 






4.4 


1.8 


LBGi 


3 


871 


9x ~ 9' X 9' 


i?. < 25.5 




7 






3.0 


1.98 ±0.3 


LBG2 


3 


268 


6x ~9 


' X9' 


R< 25.5 




7 






3.0 


1.8 



4 INTERPRETATION 

Given the general trend shown in Figure 1, we have to ad- 
dress the following questions: (1) Why does the observed 
clustering vary across datasets? (2) Is it possible to scale dif- 
ferent datasets to make a self-consistent dataset, and then 
study evolution within that superset? (3) Can we constrain 
any of the models using this data? 

There are significant discrepancies between the ampli- 
tude of clustering in different surveys even where they sam- 
ple the same redshift intervals. To some extent these discrep- 
ancies are intrinsic because different surveys sample different 
populations of galaxies. Some differences are introduced by 
the extent to which the luminosity function is probed, i.e. in 
one survey the limiting magnitude may allow one to probe 
galaxies much fainter than L* and in another case the lim- 
iting magnitude may be comparable to L*. Since different 
surveys measure clustering at different scales, and many of 
them extrapolate the values for ro{z) on scales bigger than 
the areas actually covered by the surveys themselves (see 
table 1), further differences are introduced if the galaxy cor- 
relation function is not a true power law in the range of 
scales between the scale of measurement and 8/i~^Mpc. 

All these effects are very likely to "bias" the different 
measurements with respect to each other. In this section 
we tackle these issues one by one, in order to correct the 
data for their effects so to obtain "unbiased" sub-sets of 
measurements compatible with each other. 



4.1 Scale Dependence 

As already mentioned in the previous sections and as shown 
in table 1, different surveys cover areas of the sky which vary 
greatly in size from one survey to another. This problem of 
sampling objects on different scales is likely to introduce a 
relative bias amongst clustering measurements coming from 
different surveys. In fact many of these studies quote values 
for the clustering length ro(z) which have been obtained by 
extrapolating the power-law trend of ^(r, z) to scales much 
larger than the physical scales of the surveys. In order to 
correct for this effect we write the bias as 



b (r,z) 



^m{r,z)' 



(8) 



where ^m and ^9 respectively are the mass-mass and galaxy- 
galaxy correlation function and f is some fiducial scale 
length. Note that by writing the bias as an explicit func- 
tion of f which is different for each survey, this approach 
corrects for the "scale effect", since it compares theoretical 
quantities and measurements evaluated at the same scale. 

We start by evaluating the mass-mass correlation func- 
tion £,m{r, z). Throughout this section we will use the nota- 
tions adopted in Peacock (1997). We start with a (dimen- 
sionless) primordial power-spectrum of the form Af^^ (k) oc 
fc''+^r| (with n = 1 for CDM models). The transfer func- 
tion for CDM family of models is the one given by Bardeen 
et al. (1986): 



T(fc) 



ln(l + 2.34q) 
2.34g 



[1 + 3.89g + (16.1q)^ + (5.46g)^ + (e.?!?)"*] 



-1/4 



(9) 
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where q = {k/h)/V Mpc^^ and the shape parameter F is re- 
lated to the present-day matter parameter Oq and t he bary- 
onic fraction Jlot via T = Q,q h exp [f2ob — yh/O.S r^ob/^o] 
(Sugiyama, 1995). The normahzation of the power-spectrum 
is fixed by specifying the values of erg" (see section 2), ac- 
cording to the expression: 



/ liiiN2 



AL(fc)'^^ 



9 



k {kRf 



[sin kR — kR cos kR] 



(10) 



In order to go from the linear power-spectrum Afij,(fc) to 
the non-linear case, Peacock & Dodds (1994, 1996), following 
an approach originally introduced by Hamilton et al. (1991), 
assumed there is a universal fitting function F relating the 
two according to the expressions: 



A^(fc)=F[AL(ko) 
fco= [l-l-A'(fc)] 



-1/3 



fc, 



(11) 



where fco and k are respectively the linear and non-linear 
wavenumber and F is given by: 



F{x) 



1 -f Bjlx -f [Ax] 



all 



l + (\AxYg{Q.oy^/[Vx^/^\Y 



1//3 



(12) 



(see Peacock & Dodds (1996) for the values of the param- 
eters in equation |l2[ ). g{0.o) is a suppression factor which 
measures the rate of growth of clustering in generic cosmolo- 
gies relative to the growth in an Einstein-de Sitter Universe. 
Lahav et al. (1991) and Carrol, Press & Turner (1992) found 
that this quantity can be almost exactly approximated by 



g{Q.o,Q,A,o) = -ilo 



a 



,4/7 



- f2A,0 



- f^o/2) 



(1 + 



^A.O , 

70 ' 



(13) 



The above formulae were given in Peacock & Dodds 
(1996) and Peacock (1997) for z = Q. However Moscardini 
et al. (1997) argue that they all apply to any cosmic epoch z, 
as long as one replaces g by g{z), and interpret the quantity 
X as the linear power-spectrum at epoch z: 



X = AL(fco, z) = AL(fco)(i + z)'^ [g{z)/g{Q)f 



(14) 



Replacing g by g{z) implies writing both Q. and Q.a as a 
function of redshift. Their expressions have been taken from 
Lahav et al. (1991) and obviously depend on the underlying 
cosmology: 

3 



Q 



(1 



nA(z) 



(15) 



f^o {h/ho) 

" 3/i2' 
where A is the cosmological constant which does not vary 
with time, and 

h{z) = ha [Vto{l + zf - {0.0 + Q.K,o - 1)(1 + zf + Q.k,o\ ^''' (16) 

Now that we have all the theory in place, we can finally 
evaluate the evolution of the spatial two-point correlation 
function ^(r, z') with time. Note that ^(r, z) is related to the 
power-spectrum A (fc, z) via the expression: 



has been obtained for all the redshifts sampled by differ- 
ent surveys. The fiducial length f has been chosen to be 
Tmax/'i, where rmax is the upper limit of the range of scales 
used to measure the clustering signal within each survey. 
This choice, even though somewhat "ad hoc", seems plau- 
sible given that the measurements go from about r = to 
Tmax, and this gives an effective radius ^ Vmax/'i- A more 
exact analysis would have to take into account the actual 
bin steps, weighting and fitting procedure used in deriving 
^ for each survey, but since we were only interested in the 
relative scalings, we kept it to this simple approximation. 
Note that the values of r^ax (and therefore of f) vary both 
with redshift and cosmology and correspond to scales that 
range from linear to highly non-linear (see Figure 2). In the 
cases where the 3-d clustering was obtained by deprojecting 
the angular (2-d) correlation function w{9) (as for instance 
in the HDF points), the value of rmax was obtained via 
r — 9maxX, with Omax maximum angular scale and x comov- 
ing coordinate, whose expressions for different cosmologies 
are given in Appendix A. We repeated this analysis for the 
three cosmologies introduced in section 2. 

Figure 2 shows the theoretical predictions for the spa- 
tial two-point correlation function ^(f, z) at different red- 
shifts and compares them with the respective data. Each 
panel corresponds to a different survey (an exception has 
been made for the HDF in order not to overcrowd the rel- 
ative panel) and the different curves show the predictions 
for ^ at the redshifts z quoted for each of their clustering 
measurements; the lower curves are for higher z's. 

The points show ^(f , z) derived from the observations 
assuming a power-law form ^{f,z) = [f/ro(z)]~^. The val- 
ues of 7 are directly observed, as listed in table 1 and ro{z) 
is derived for each cosmology from expressions (0) and (m. 
Note that the assumption of a power-law form for ^ is justi- 
fied by the small areas covered by most of the surveys. In the 
case of wide-area surveys (such as Stromlo-APM and IRAS) 
f has been fixed to a standard value of 10 Mpc, indepen- 
dent of cosmology. Different symbols correspond to different 
surveys. In the HDF panels crosses indicate the results from 
Magliocchetti & Maddox (1999) (HDF2), while empty cir- 
cles are for Connolly, Szalay & Brummer (1998) (HDFi). In 
the LEG panel the higher point is taken from Adelberger 
et al. (1998). Figure 2 shows the great spread in the scales 
sampled by different surveys ranging from a few tenths to 
>10 Mpc. It is also possible to note that even within the 
same survey (e.g. HDF and LBG), different analysing tech- 
niques measure different scales. As we will show later in this 
section, this scale effect can partially explain the apparent 
discrepancy amongst results quoted for similar sets of data. 



C(r,z)= / A\k,z) 



, sin kr dk 



kr k 
In order to compare theoretical models with data, ^(f , z) 



4.2 Populations 

Different surveys in general sample different populations of 
objects; this is due to both the selection criteria and the 
redshift range sampled by the survey. Objects selected in 
the UV band will be dominated by star-forming galaxies, 
while for instance, B band selected objects will give mix of 
early and late-type galaxies. At higher redshifts the popu- 
lation mix sampled by a survey will also depend on z since 
the rest-frame pass band is shifted towards the blue. For 
example a survey which selects objects in the I-band will 
contain many early- type galaxies at low z, but for 2 > 1.5 
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Figure 2. Theoretical predictions for the spatial two-point correlation function §(-R) at different redshifts and for an Einstein-de Sitter 
Universe. Each panel corresponds to a different survey and the dashed lines indicate the models obtained for those redshifts z's which 
have available measurements of ro(S) (or alternatively of a^{z), see table 1). Lower lines correspond higher redshifts. The points show the 
values for i{K) obtained from the data in the assumption of a power-law form for the correlation function, where R is half the size of the 
maximum scale used to derive the clustering measurements, representing the scale used to determine rp. Different symbols correspond 
to different surveys, fn the HDF panels crosses indicate the results from Magliocchetti &; Maddox (1999) (HDF2), while empty circles 
are for Connolly, Szalay & Brummer (1998) (HDFi). In the LBG panel the higher point is taken from Adelberger et al. (1998). 



the observed I-band corresponds to the rest-frame UV band, 
and so the sample will be dominated by star-forming galax- 
ies. This also implies that, even within the same survey, the 
galaxy population that is sampled depends strongly on the 
redshift. 

Since late-type/star-forming galaxies cluster more 
weakly than early-type ones (e.g Loveday et al. 1995; Her- 
mit et al. 1996; Guzzo et al. 1997), the changes in the galaxy 
population must be taken into account when comparing the 
clustering measurements from different surveys and different 
redshifts. 

As a simple first step to minimizing this effect we di- 



vided our sample into different populations, according to the 
rest-frame selection band. The Stromlo-APM survey sam- 
ples galaxies in the rest B-band (see table 1) as it does the 
Las Campanas redshift survey (LCRS, objects selected in 
the observed R-band at a median redshift z — 0.13) and 
the CN0C2. The same holds for both the CFRS survey (I- 
band selected objects at 0.3 ^ z < 1) and the HDF for 
2; < 1.4. We will denote the objects selected by these sur- 
veys with the collective name of blue. At higher redshift the 
observed Lband shifts into the rest-frame UV band so that, 
for z > 1.4, the population sampled by the HDF will be 
dominated by star- forming galaxies. The selection criteria 
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of both IRAS and Lyman-Break galaxies also lead to sam- 
ples dominated by galaxies undergoing star-formation. We 
will call these objects star-forming. The Keck K-band sur- 
vey selects objects in the near-IR, which for 0.3 ^ z <^ 1.4 
corresponds to a rest-frame R-band. We therefore expect the 
populations sampled to be mainly early-type (red) galaxies. 

Now that we have made this division into different pop- 
ulations according to their star-forming activity and there- 
fore their colours, we can go back to the issues of the evolu- 
tion of galaxy clustering and the redshift evolution of bias. 
In more detail, for each of the three populations we have 
then evaluated the quantity b{f, z) = \/[Cs(^j z)/£^rn{r, z)], 
as explained in section 4.1. 

Figure 3 shows the trend of b{z) as a function of (l + z), 
respectively for blue (left panels), red (middle panels) and 
star- forming galaxies (right panels). Different symbols cor- 
respond to different surveys and are the same as in figure 2. 
Once again the upper LBG point (empty triangle) is taken 
from Adelberger et al. (1998). The upper panels are for an 
Einstein-de Sitter Universe (f2o = i, ho — 1), the middle 
panels for open cosmologies (fio = 0.4, ho = 0.65) and the 
lower panels correspond to flat geometries with non-null cos- 
mological constant {flo = 0.4, r2A,o = 0.6, ho = 0.65). 

There are many interesting features to note in Figure 
3. The first is the trend of points at low redshift {z < 1.5). 
The large spread amongst the different surveys seen in Fig- 
ure 1 is much reduced after renormalizing the amplitude 
according to the scale sampled by each survey. Panels on 
the left-hand side of Figure 3 show the corrected results in 
excellent agreement with each other. This good agreement 
also suggests that there is little dependence of bias on scale. 

The systematic higher clustering amplitude in the Keck 
K-band survey as compared to other surveys (see middle 
panels in Figure 3) can now be explained as due to different 
clustering properties of different populations. In more detail 
the results show, as expected from local measurements, that 
early-type objects cluster more strongly than late- type ones. 

In Figure 3 we have also plotted some theoretical pre- 
dictions for the evolution of bias with redshift. The solid 
lines correspond to the test-particle model (Bl), described 
by equation (0). Note that the value of the quantity bo = 
<^s,g/o's,m (where in this case as,m is the non-linear r.m.s. 
mass density fluctuation at 8h~^ Mpc, as calculated from the 
models described in section 4.2) varies for different popula- 
tions. For the models we use a different value of as,g for each 
galaxy type, with star-bursting galaxies being less clustered 
than quiescent galaxies with old stellar populations. We es- 
timate the amplitudes ag.g from the values of ro measured 
for three subsamples from Stromlo-APM survey (Loveday, 
Tresse & Maddox, 1999); for galaxies with no emission lines 
(red objects) as,g ~ 1.13; galaxies with weak emission (blue 
galaxies) as,g = 0.93; and as,g = 0.66 for galaxies with 
strong emission lines (star- forming galaxies). 

The dashed lines correspond to the merging model for 
the evolution of bias with z (B2), with lower curves corre- 
sponding to lower halo-masses (IO^^Mq < Mmin < W^'^Mq, 
see section 2). Note that neither the Bl nor the B2 curves 
have been anchored to any z — point. 

At low redshift, quite independent of the cosmology, 
the Bl models seem to provide a reasonably good fit to 
the data, both in the case of blue and early-type objects, 
even though we cannot rule out the constant bias model 



as a possible description (BO - equation hi, lines parallel to 
the X axis, not shown in the figure). This is also seen in the 
analysis of the clustering of radio galaxies (e.g. Magliocchetti 
et al. 1999); in this case the mean redshift for the clustering 
measurements is about z ~ 1. Even though the data cannot 
distinguish between constant bias (BO) models and models 
with bias linearly evolving with redshift (Bl), it is clear that 
the clustering measurements obtained for the radio sample 
in such a redshift range are in confiict with the predictions 
obtained from the merging (B2) model. 

At higher redshifts {z > 2) though, both the constant 
bias and the test particle model grossly fail to describe the 
data which show very high values for b{z). Note that in this 
redshift range we have measurements only for star-forming 
objects, therefore it is impossible to state whether this rise 
in the level of biasing is always true for z > 2, independent of 
the population, or only holds for star-forming galaxies. Al- 
though the match is not perfect, the merging model correctly 
describes the trend of the data (and also matches the IRAS 
points), with lower halo-masses {Mmin = 10^^ Mq) required 
for an Einstein-de Sitter Universe and higher halo-masses 
(Mmin ~ W^^Mq) for low-deuslty models (independent of 
f^A.o)- 

4.3 Malmquist bias 

Now we turn to a discussion of the selection effects, namely 
the Malmquist bias and K-correction. 

Although it is not possible to correct for the Malmquist- 
bias and K-correction, models can be made to allow for 
their effects. In a detailed model (e.g. by using the semi- 



analytic approach; (Kauffmann et al. 
I99S|; 



Somerville fc Primack 1998 )), these selection effects 



1999; Baugh et al 



can be convolved with the model and direct comparison can 
be made with the observations. However, these models con- 
tain a large number of parameters and the gain in details 
is accompanied by a loss of intuitive understanding. In or- 
der to keep the discussion simple and independent of model 
parameters, we will not use this approach here. 

Instead, by dividing the objects into different popula- 
tions we minimize the K-correction problem, because objects 
belonging to the same population are likely to have similar 
spectra, and so have similar K-corrections. 

As already mentioned in the introduction, more lumi- 



nous galaxies found to clus ter more strongly (Park ot al 
1994; Loveday et al. 1995) than fainter ones. Therefore, 
for the apparent magnitude selected samples in Figure 3 
we would expect to see an increase of the clustering in- 
troduced by the Malmquist bias at higher z's. The CFRS 
measurements are based on all galaxies brighter than the 
apparent magnitude limit of the survey, and so the points 
sample objects with increasing absolute magnitudes as one 
goes to higher redshifts. We can see a hint of the expected 
Malmquist effect in these points in that the third bin (cen- 
tred at z ~ 0.87) shows a slight increase in clustering ampli- 
tude, and includes only galaxies brighter than L*g (see table 
1). However this is at a very low significance level. On the 
other hand the measurements from the CNOC2 and the Las 
Campanas surveys are based on sub-samples of galaxies with 
a cut-off in absolute magnitude (respectively Lr > —20 and 
—23 < Lr < —18, and so should not show this effect. 

In the Keck K-band survey (panels in the centre of Fig- 
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Figure 3. Trend of the bias 6(2:) as a function of redshift for different populations of objects. Left panels correspond to blue objects, 
middle panels to red objects and right panels to star-forming galaxies. Symbols are the same as in figure 2. Results have been obtained 
for an Einstein-de Sitter Universe (top), open Universe (centre) and flat Universe with a cosmological constant (bottom). Lines indicate 
theoretical models for the evolution of bias as explained in section 2. Solid lines represent the test-particle (Bl) model, while the dashed 
lines are for the merging (B2) model with Mmin = 10^^, 10^^ and IO^^Mq, lower curves corresponding to lower halo-masses (see text 
for details). 



ure 3), clustering measurements have been obtained for all 
the galaxies respectively with Lk > —23 in the three higher 
redshift bins and with Lk > — 2L5 in the lower redshift 
range. We therefore expect the lower-2 point to underesti- 
mate the clustering with respect to the results for higher 
z's and brighter absolute magnitudes. However, given that 
Lk "^ ^25 + 5 log(h5o), even for Lk > —23.5, all the galax- 
ies still sit on the flat faint-end of the luminosity function. 
Since Malmquist bias is likely to be important only for lumi- 
nosities > L*, we argue that the effect should small in this 
case. 

Unfortunately we do not have reliable information on 
the absolute luminosities of objects in the HDF because the 



K-corrections are highly uncertain, particularly for 2 > 1. 
However, at low redshift {z < 1.4), the faint apparent mag- 
nitude limit suggests that all galaxies will be below L* , and 
so Malmquist bias should not play an important role. This 
is consistent with the trend shown by the data in Figure 3 
(left panels). At higher redshift the galaxies in the sample 
are brighter than L*, therefore we might expect Malmquist 
bias to be more important, so that the clustering amplitude 
would increase with increasing redshift as seen in the data. 
It would be possible to quantify the Malmquist bias ef- 
fect by converting the absolute luminosity limit for each sub- 
set of galaxies to an equivalent mass limit and hence find the 
expected clustering amplitude for each redshift. However the 
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Figure 4. The estimated comoving galaxy space density as a 
function of redshift. The points show the estimates for the HDF 
galaxies as used in the clustering measurements. The hnes show 
the Press-Schechter prediction for the density of halos with mass 
greater than Af„i„ = 10^°, lO", lO^^ and lO^^jv/g, as labelled. 
We applied the standard Press-Schechter formalism using the cos- 
mology and power spectrum from our model I. The data at z > 3 
are almost certainly affected by random errors, because the space 
density implies a low M^i„ , whereas the clustering implies a very 
high Mmin (see Figure |). 



K- correct ions, and hence the derived absolute magnitudes, 
are highly uncertain at these redshifts. Also the mass-to- 
light ratio is unknown at high redshift, and given the rapid 
evolution of stellar populations at z < 1, it is likely to be 
very different to local galaxies. An alternative way to tackle 
this problem is to consider the space density of objects. In 
models based on high peak biasing, rarer objects correspond 
to higher masses, which are more strongly clustered. Hence 
a comparison between the observed and predicted comoving 
space-density of galaxies allows us to estimate the effective 
Mmin for each redshift bin, independent of uncertainties in 
the absolute magnitudes. This should provide a more robust 
approach to quantifying the Malmquist bias. 

In Figure m the points show the observed mean density 
for each HDF redshift bin, estimated as iV = Naz/V, where 
V is the comoving volume in each bin assuming an Einstein- 
de Sitter Universe. Note that, as expected, there is a clear 
trend for the space-density to decrease with increasing red- 
shift, although the decrease is not very strong, even at high 
redshift. The lines show the space density of haloes of mass 
greater than Mmin = 10^°, lO", 10^^ and IO^Mq predicted 
from the standard Press-Schechter formalism with the power 
spectrum as in our model I. For z ^ 2 the points are con- 
sistent with a fixed halo mass ~ 5 x 10^° Mq. For 2 < z < 4 
there is a hint that Mmin rnay rise to ~ 10^^ Mq, but the 
uncertainties are large, and for 2 > 4 the apparent decrease 
is not significant. The merging model (B2) introduced by 
Matarrese et al. (1997) and Moscardini et al. (1998) predicts 
the evolution of bias with redshift, according to the mass of 



the halos. Thus we can compare directly to the model (B2) 
(dashed lines in the panels on the right-hand side of Figure 
3) with a minimum halo mass Mminiz) which is roughly 
constant or slightly increasing with at higher redshifts (see 
also Arnouts et al. 1999). 



5 CONSTRAINTS ON «„ 

The rate of evolution of the galaxy correlation function with 
r edshift has be en used to estimate Cosmological parameters 
(Peacock 1997). This approach assumes that the bias does 
not have any scale dependence, and, it does not vary sig- 
nificantly in time. However, recent studies have questioned 
this approach and have shown that the time evolution of 
bias is very important, even if we can choose to ignore its 
scale dependence at large scales. The most remarkable man- 
ifestation of the time evolution of bias are the Lyman break 
galaxies, that have a clustering amplitude comparable to 
present day galaxies. The models of bias evolution discussed 
in section 2 also suggest that the time evolution of bias is 
very important. 

In principle, it is possible to use the models of bias along 
with the observations to see if the observed redshift evolu- 
tion of clustering is consistent with a given model or not. 
However, this approach requires the power spectrum for den- 
sity fluctuations as an input, and hence is model dependent. 

We propose a simple test to derive a lower limit on flo- 
We make use of the fact that in all realistic scenarios, bias 
always increases as we go to higher redshifts. This monotonic 
increase in bias implies that the galaxy correlation function 
will evolve at a slower rate than the correlation function of 
dark matter, i.e. Dgai{z) < Dm{z). Dm{z) depends only on 
Jlo (and JIa.o), and evolves at a slower rate for lower flo- If 
we compute the allowed values of Dgai{z) in a given survey, 
and we get a firm lower limit on this rate, then we can rule 
out those cosmologies that predict Dm{z) < DgaLobsv{z). 
We would like to stress, that this method is completely model 
independent if used at large/linear scales. 

To further illustrate this idea, we have plotted the best 
fit points for CFRS and CNOC2 surveys on the flo-os plane 
in figure U. Thick lines show the confidence limits at the 
68% level. The error bars at this confidence level are too 
big to be of much use. We have not shown data from any 
other survey listed in this paper because the error bars for 
those are even larger. However, to demonstrate what can be 
achieved by the next generation surveys, like 2DF (Maddox 
1998), Sloan (Kim et al. 1999) and VIRMOS (Garilli et al. 
1999) which will have many more galaxy redshifts, we have 
shown the contours for error bars reduced by a factor 5. In 
this case, one may begin to rule out interesting regions of 
the parameter space. 

This test is more effective at lower redshifts as the dif- 
ference in the growth rate in different cosmologies is more 
striking. Thus surveys like Sloan and 2DF may provide some 
useful constraints through this test. 
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Figure 5. This figure shows the best fit values of Hq and trg g(0) 
for CFRS and CNOC2 surveys assuming an open universe with 
f^A 0=0- The 68% confidence limits are shown as thick lines. 
Dashed lines show the same if error bars are reduced to 1/5 of 
their present size. Allowed values for Qq lie on the right-hand side 
of the diagram. See text for details. 



6 CONCLUSIONS 

We have presented a detailed analysis of the evolution of 
galaxy clustering with redshift, specifically aimed at con- 
straining epoch-dependent models for the bias. 

As a first step we put together all the available measure- 
ments from different surveys, covering z ~ to z ~ 4.5, and 
we calculated the values for the r.m.s galaxy fluctuation den- 
sity on a 8h~^ Mpc scale, o-g.g. This has been repeated for 
different cosmologies. We find that the amplitude of galaxy 
clustering shows a general trend of decreasing at redshifts 
z < 2 and increasing beyond that. 

The amplitudes for different surveys show a remark- 
able scatter between each other, and there are several fac- 
tors which might be responsible for this. First is the "scale- 
dependence" problem, coming from the fact that the var- 
ious surveys measure clustering at different scales. Sec- 
ond is that different surveys select different populations of 
objects which do not have similar clustering amplitudes. 
Third is the Malmquist bias, since brighter objects clus- 
ter more strongly than fainter ones. By using the Peacock 
& Dodds (1994, 1996) models for non-linear evolution of 
the power-spectrum P{k, z) (and consequently of the two- 
point correlation function ^(r, z)), we corrected for the scale- 
dependence effect and obtained values for the bias written as 



h\f, 



■ ^g{r, z)/^rnir, z), where r is half of the maximum 



scale used within each survey to work out estimates on the 
clustering amplitude. Note that this approach can in princi- 
ple also show a dependence of bias on scale. We allowed for 
the population selection effect by dividing the objects into 
three categories, according to the rest-frame pass band of 
their selection (roughly corresponding to their star-forming 
activity) and analyzing the evolution of clustering for each 



of these populations separately. We also considered the ex- 
pected variations introduced by the Malmquist bias. 

We find that the main cause for the scatter amongst 
different surveys is the scale-dependence problem. Renor- 
malizing the amplitudes according to the physical scale at 
which each measurement was taken brings the results in ex- 
cellent agreement with each other. 

The main conclusion of this work is that the bias grows 
monotonically from the present epoch to high redshifts, and 
the rate at which bias grows increases rapidly as we go to- 
wards higher redshifts. 

At low redshift the trend of the data-points seems to 
suggest a bias whose functional form is a weakly increasing 
function of the redshift. The test-particle bias is a useful 
model that allows us to predict the evolution of bias for 
objects that are positively biased. We do not see strong evi- 
dence for scale dependent bias or Malmquist effects, but the 
uncertainties are large. 

At higher redshifts the clustering signal appears to come 
from objects which are highly biased with respect to the un- 
derlying distribution of mass. The merging model for the 
evolution of bias (B2) can correctly describe the steep rise 
in clustering amplitude seen for 2 > 2. We note however 
that there are several uncertainties in the interpretation of 
clustering at such redshifts. In particular, at high z we have 
measurements only for the population of star-forming galax- 
ies, and the lack of reliable absolute magnitude estimates for 
these objects makes it impossible to apply any quantitative 
correction for the Malmquist bias. 

Lastly we find that the rate at which the observed am- 
plitude decreases at low redshifts is slower than the linear 
rate of evolution for density perturbations in dark matter in 
most models. We present arguments which allow us to rule 
out all models for which this is not the case, i.e. the rate of 
evolution of galaxy clustering is faster than the linear rate. 
Present observations do not rule out any model but future 
observations will allow us to constrain the density param- 
eter. We would like to emphasise that this test does not 
depend on any detailed modelling of galaxies/halos and so 
should provide a reliable constraint. 
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APPENDIX A: USEFUL COSMOLOGICAL 
QUANTITIES 

We express here the quantities appearing in equation 0. Note 
that P{Q,o, z) is obtained by direct integration of the Limber 
equation (Peebles, 1980). 
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